Search for: All records

We use numerical simulations of circumplanetary disks to determine the boundary between disks that are radially truncated by the tidal potential and those where gas escapes the Hill sphere. We consider a model problem, in which a coplanar circumplanetary disk is resupplied with gas at an injection radius smaller than the Hill radius. We evolve the disk using thePhantomsmoothed particle hydrodynamics code until a steady state is reached. We find that the most significant dependence of the truncation boundary is on the disk aspect ratioH/R. Circumplanetary disks are efficiently truncated forH/R≲ 0.2. ForH/R≃ 0.3, up to about half of the injected mass, depending on the injection radius, flows outward through the decretion disk and escapes. As expected from analytic arguments, the conditions (H/Rand Shakura–Sunyaevα) required for tidal truncation are independent of planet mass. A simulation with largerα= 0.1 shows stronger outflow than one withα= 0.01, but the dependence on transport efficiency is less important than variations ofH/R. Our results suggest two distinct classes of circumplanetary disks: tidally truncated thin disks with dust-poor outer regions, and thicker actively decreting disks with enhanced dust-to-gas ratios. Applying our results to the PDS 70 c system, we predict a largely truncated circumplanetary disk, but it is possible that enough mass escapes to support an outward flow of dust that could explain the observed disk size.

 

Young protostellar discs are likely to be both self-gravitating, and to support grain growth to sizes where the particles decoupled from the gas. This combination could lead to short-wavelength fragmentation of the solid component in otherwise non-fragmenting gas discs, forming Earth-mass solid cores during the Class 0/I stages of young stellar object evolution. We use three-dimensional smoothed particle hydrodynamics simulations of two-fluid discs, in the regime where the Stokes number of the particles St > 1, to study how the formation of solid clumps depends on the disc-to-star mass ratio, the strength of gravitational instability, and the Stokes number. Gravitational instability of the simulated discs is sustained by local cooling. We find that the ability of the spiral structures to concentrate solids increases with the cooling time and decreases with the Stokes number, while the relative dynamical temperature between gas and dust of the particles decreases with the cooling time and the disc-to-star mass ratio and increases with the Stokes number. Dust collapse occurs in a subset of high disc mass simulations, yielding clumps whose mass is close to linear theory estimates, namely 1–10 M⊕. Our results suggest that if planet formation occurs via this mechanism, the best conditions correspond to near the end of the self-gravitating phase, when the cooling time is long and the Stokes number close to unity.

 

Diversity in the properties of exoplanetary systems arises, in part, from dynamical evolution that occurs after planet formation. We use numerical integrations to explore the relative role of secular and resonant dynamics in the long-term evolution of model planetary systems, made up of three equal mass giant planets on initially eccentric orbits. The range of separations studied is dominated by secular processes, but intersects chains of high-order mean-motion resonances. Over time-scales of 108 orbits, the secular evolution of the simulated systems is predominantly regular. High-order resonant chains, however, can be a significant source of angular momentum deficit (AMD), leading to instability. Using a time series analysis based on a Hilbert transform, we associate instability with broad islands of chaotic evolution. Previous work has suggested that first-order resonances could modify the AMD of nominally secular systems and facilitate secular chaos. We find that higher order resonances, when present in chains, can have similar impacts.

 

ABSTRACT The dissociation and ionization of hydrogen, during the formation of giant planets via core accretion, reduce the effective adiabatic index γ of the gas and could trigger dynamical instability. We generalize the analysis of Chandrasekhar, who determined that the threshold for instability of a self-gravitating hydrostatic body lies at γ = 4/3, to account for the presence of a planetary core, which we model as an incompressible fluid. We show that the dominant effect of the core is to stabilize the envelope to radial perturbations, in some cases completely (i.e. for all γ > 1). When instability is possible, unstable planetary configurations occupy a strip of γ values whose upper boundary falls below γ = 4/3. Fiducial evolutionary tracks of giant planets forming through core accretion appear unlikely to cross the dynamical instability strip that we define. 

ABSTRACT The core accretion model of giant planet formation has been challenged by the discovery of recycling flows between the planetary envelope and the disc that can slow or stall envelope accretion. We carry out 3D radiation hydrodynamic simulations with an updated opacity compilation to model the proto-Jupiter’s envelope. To isolate the 3D effects of convection and recycling, we simulate both isolated spherical envelopes and envelopes embedded in discs. The envelopes are heated at given rates to achieve steady states, enabling comparisons with 1D models. We vary envelope properties to obtain both radiative and convective solutions. Using a passive scalar, we observe significant mass recycling on the orbital time-scale. For a radiative envelope, recycling can only penetrate from the disc surface until ∼0.1–0.2 planetary Hill radii, while for a convective envelope, the convective motion can ‘dredge up’ the deeper part of the envelope so that the entire convective envelope is recycled efficiently. This recycling, however, has only limited effects on the envelopes’ thermal structure. The radiative envelope embedded in the disc has identical structure as the isolated envelope. The convective envelope has a slightly higher density when it is embedded in the disc. We introduce a modified 1D approach which can fully reproduce our 3D simulations. With our updated opacity and 1D model, we recompute Jupiter’s envelope accretion with a 10 M⊕ core, and the time-scale to runaway accretion is shorter than the disc lifetime as in prior studies. Finally, we discuss the implications of the efficient recycling on the observed chemical abundances of the planetary atmosphere (especially for super-Earths and mini-Neptunes). 

ABSTRACT In dense star clusters, such as globular and open clusters, dynamical interactions between stars and black holes (BHs) can be extremely frequent, leading to various astrophysical transients. Close encounters between a star and a stellar mass BH make it possible for the star to be tidally disrupted by the BH. Due to the relative low mass of the BH and the small cross-section of the tidal disruption event (TDE) for cases with high penetration, disruptions caused by close encounters are usually partial disruptions. The existence of the remnant stellar core and its non-negligible mass compared to the stellar mass BH alters the accretion process significantly. We study this problem with SPH simulations using the code Phantom, with the inclusion of radiation pressure, which is important for small mass BHs. Additionally, we develop a new, more general method of computing the fallback rate which does not rely on any approximation. Our study shows that the powerlaw slope of the fallback rate has a strong dependence on the mass of the BH in the stellar mass BH regime. Furthermore, in this regime, self-gravity of the fallback stream and local instabilities become more significant, and cause the disrupted material to collapse into small clumps before returning to the BH. This results in an abrupt increase of the fallback rate, which can significantly deviate from a powerlaw. Our results will help in the identification of TDEs by stellar mass BHs in dense clusters.